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PERIODIC AND ASYMPTOTIC ORBITS FIVE BODY 


DANIEL BUCHANAN? 


Abstract 


Periodic and asymptotic orbits are obtained for the motion infinitesimal 
body subject the Newtonian attraction four finite bodies. The finite 
bodies are equal mass and remain relatively fixed the vertices square 
while they revolve about their common centre gravity with uniform angular 
velocity. The infinitesimal body moves the vicinity the libration 
points the finite bodies. Asymptotic orbits two dimensions and periodic 
orbits two and three dimensions are obtained. 


Introduction 


This paper deals with the determination periodic and asymptotic orbits 
the vicinity the libration points restricted problem five bodies. 
Four the bodies are equal finite masses and are assumed fixed 
the vertices square that rotates uniformly about the centre gravity 
the four bodies. The fifth body attracted the finite 
bodies according the Newtonian law the inverse square. Periodic orbits 
both in, and cutting, the plane motion the finite bodies and asymptotic 
orbits within the above plane motion are obtained for the infinitesimal 
body moving the vicinity its points libration. 


The Libration Points 


recent paper Hinrichsen (17) considered the libration equilibrium 
points the general problem bodies which one body infinitesimal 
and the other bodies are equal finite masses. The finite bodies are 
relatively fixed the vertices regular 1)-sided polygon and are 
rotating uniformly about their common centre gravity. 

The following axes and units were selected Hinrichsen. The motion 
the infinitesimal body referred set rectangular co-ordinate axes 
that rotate uniformly with the finite bodies and have the origin their 
common centre gravity. The x-axis chosen that, for even, one 
vertex the polygon lies the negative half the x-axis, and for odd 
the vertices lie both positive and negative sides the x-axis. The co- 
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and those the infinitesimal (x, The unit mass chosen the 
sum the masses the finite bodies and the unit distance the distance 
from any the finite bodies the centre gravity the other bodies. 
The unit time taken that the angular velocity unity. 


With the axes and units thus chosen, the differential equations motion 
for the infinitesimal body are 


(1) 
“dt Q,, 
Q 1 B +2) 
+2, 


These equations admit the Jacobi integral 


The libration points, found Hinrichsen, following the methods 
Lagrange (18, Vol. pp. 229-324; 20, pp. 309-318) and others, are the 
solutions 


n-1 1 
(1 - 4) @ - 209 = Q, 

n-1 (2) 


multiplying the first and second these equations y,) and 
respectively, and subtracting obtain 


These equations denote the various lines through the vertices the 1)- 
sided polygon and also the circles unit radius having their centres the 
vertices the polygon. These lines and circles divide the plane into various 
For the largest value considered detail Hinrichsen, the plane 
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divided Fig. the libration points being the shaded regions. They 
are found follows: 


(0, 


all there are points, all points being within distance 7/4 from the 
origin. 
The Results Obtained 
Owing the symmetry the problem, orbits only the vicinity the 
following are considered: 


The origin, 


(a) (1.255, 0), 
(b) (0.465, 0.465), 
(c) (0.625, 0.625). 


The corresponding orbits near the nine other points can obtained inter- 
changing changing the signs and and the co-ordinates (a), 

The method used finding the periodic orbits follows somewhat closely 
that used Moulton (21, Chap. his treatment analogous three 
body problem. use, however, the operator differential equations 
instead the more laborious transformation normal form used 
Moulton. 

Similar methods were used Buck (16, Chap. IX), Buchanan (14, 15), 
and the author (1, Chap. 10, 2-6)* treating various problems periodic 
orbits. 


*Other articles, chiefly the the Royal Society Canada. 
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The asymptotic orbits are found the method used the present author 
certain other analogous problems (7-13)*. 

Periodic orbits two dimensions are obtained for motion near (a) and 
These orbits correspond orbits Class (21, Chap. V). 
Asymptotic orbits two dimensions are found near the four points, the 
origin, (a), (6), and (c). these orbits the infinitesimal body approaches 
the equilibrium points the time approaches infinity. 

For motion outside the xy-plane, periodic orbits for all four points are 
obtained. The orbit near the origin, however, linear, being along the 
z-axis. These orbits correspond Moulton’s orbits Class 

Orbits category Class were not considered. They depend 
upon the commensurability two parameters, and their determination even 
when this commensurability assumed quite complex. 

Asymptotic orbits three dimensions similar those obtained the 
author other cases are deferred for future treatment. 


The Equations Displacement 


order discuss the motion the vicinity the equilibrium points 
substitute Equations (1) 
where (A, are the co-ordinates equilibrium point; are new depend- 
ent variables; the new dependent variable; parameter denoting the 
scale factor the orbits; and parameter vanishing with sub- 
stituting Equations (5) Equations (1) obtain 


/ 


The asterisk denotes the omission terms independent ev. 
consider only the terms Equations (6) that are linear and 
obtain, after dropping the equations variation, 


(7) 


*Other ar.icles, chiefly the Transactions the Roval Society Canada. 
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The values for the and d’s for the equilibrium points considered are listed 


Table 
TABLE 
point (0, (0.465, 0.465) (0.625, 0.625) 
| 3 | 
Adé 3 | | 
A} 3 | 
A; 3 | | 
2.185 5.140 2.654 
The indicial equation Equations (7) 
Its roots for the various equilibrium points are listed Table II. 
TABLE 
Libration point Indicial equation Roots 
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Now the oscillations near libration point are stable and periodic for those 
roots the corresponding indicial equations that are purely imaginary, 
and unstable and asymptotic for the roots that are real complex. Hence 
there are asymptotic orbits near all the libration points but periodic orbits 
near (1.255, and (0.465, 0.465) only. 


Periodic Orbits the Plane 


The periodic orbits will constructed the construction similar 
for both points (1.255, and (0.465, 0.465), only the latter point will 
considered detail. 


this case the solutions the equations variation are 
(—0.61 0.79 sin Br, 


— 


: 2a d. 2a 


The trigonometric rather than the exponential form used for the periodic 
solutions order simplify the construction. 


For sufficiently small numerical values the terms Equations (6) may 
expanded converging series. The expansion yields, after has been 
cancelled off, 


x 
J 


tel 


far the computation has been carried out have 


wm 


yl 


mle 


0 
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4 1 

4 1 5 


For the equilibrium points having 
(1) (1) 1) 1 


then 

Now let 
7=0 7=0 


substituted Equations (10) and let the resulting equations cited 
Equations equating the various powers Equations 
obtain sets differential equations defining the and v;. will 
shown that the and the arbitrary constants each step the integration 
can determined that and will periodic with the period 27/8 and 
will satisfy certain assigned initial conditions. Because symmetry will 
unity inasmuch and carry the arbitrary factor Equation With 


The differential equations this step are the same the equations 


variation. Their solutions are therefore the same Equations (9). When 
the periodicity and initial conditions are imposed obtain 


(13) 
cos 0.34 sin Br, 
the second step the differential equations are 


2.99 4.24 cos 3.53 sin 287, 


| 
= 
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The complementary functions having the period 27/8 are 
where and are constants integration. 
The particular integrals, symbolic form, are 


= A | 


Since vanishes for the terms cos and sin the numerators 
Equations (16) must vanish. These terms are 


and (17) 


The determinant the terms the brackets vanishes for and 
hence one equation Equations (17) consequence the other. Selecting 
the first equation must have vanishing coefficients for both cos and 
sin Br. 
Thus 


Now the determinant 3.50 and —1.45 these two equations 
vanishes and hence both equations are satisfied when either equation is. 
From the first equation obtain With this value the parti- 
cular integrals become 


\ 


2.42 0.09 cos 287 0.06 sin 287, 
2.42 0.09 cos 287 0.06 sin 


Combining these particular integrals with Complementary Functions (15) 
obtain the desired solutions 
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2.42 (2.33 cos 0.80 sin Br) 
—0.09 cos 0.06 sin 287, 


2.42 (2.33 cos 0.80 sin Br) 


—0.09 cos 287 0.06 sin 


necessary consider detail one more step the integration before 
proceeding the general The differential equations the third step 
are 


(20) 


(3) + (3) rr . . 
where are known constants. The terms containing factor 


occur only written. order that the particular integrals 
(20) may periodic must have, Equations (17), 


the previous step the functional determinant the terms the brackets 
vanishes for Hence one equation Equations (21) con- 
sequence the other. Selecting the first equation, must have vanishing 
coefficients for both cos and sin Br. Thus 


The determinant the terms the brackets vanishes for 1.50 
and hence one equation consequence the other. Solving the first one 


for obtain 


26.45 


| 
| 
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With thus determined the particular integrals will have the same form 
and the terms cos and sin are omitted. Thus the 
complete solutions which are periodic and satisfy the initial conditions have 


the form 
3 . . . 


= u»(—T), 


(3 


3 
K 3 


7=0 
The General Step 

The remaining steps the integration are similar the preceding step. 
solving for and will necessary and possible determine 
that the coefficients cos and sin the equations analogous 
Equations (21) will vanish. The coefficients these equations are the 
same those the third step, and hence can chosen that the 
solutions will periodic. When the initial and periodicity conditions are 


imposed, the solutions take the form 


n+1 . . . 
n+l 


substituting Equations (23) Equations (11) and the results Equations 
(5) the final form the periodic solutions near 


Approximate Periodic Orbit 
consider the linear terms Periodic Solutions (24) obtain 


This represents ellipse with the centre the equilibrium point and the 
axes rotated through The semi-axes are and 0.34 
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initial values other than v(0) were chosen, say u(0) 
v(0) the approximate orbit would still ellipse oriented Fig. 
the difference being that the initial projection the latter case would not 


from apse. 
Ya 


Fic. 


Asymptotic Orbits Near (0.465, 0.465) 


construct the asymptotic orbits begin with Equations (5) but 
not change the independent variable The equations for the displacement 
from the equilibrium point (A, are the same Equations (6) except that 
neglected and considered denoting d/dt. The generating solutions 
are, Equations (9), 


where the periodic terms have been dropped. 

Now substitute Equations (11), with Equations (10) and denote 
the resulting equations Equations The solutions for are then 
found integrating Equations step step and retaining only the 
terms powers choose the initial conditions u(0) The 


desired solutions the first step are therefore 


-at 


The differential equations the second step are 


where and are known constants. The complete solutions these 
equations that satisfy the asymptotic and initial conditions are 

= 4. ea’, 


(2 (2; (2) (2) 


/ 
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The remaining steps the integration are similar this step. there 
are terms the right members the differential equations any 
step, there will terms the particular integrals. The solutions 


the general step have the form 


(i+1) -kat 


substituting Equations (25) Equations (11) and the results Equations 
(15) obtain 


k=1 


Considering only the linear terms find approximate orbit 


t 


This represents straight line through the equilibrium point having the slope 


Orbits Near (1.255, 


Only the results for the orbits near the equilibrium point (1.255, will 
listed. 


The Periodic Orbit 
The initial conditions are chosen that the infinitesimal initially projected 
from the x-axis and distance from the equilibrium point. The solutions 


are 
7 
j7=2 k= 


where are constants. The approximate shape the orbit found 


‘in Fig. ellipse with the centre the equilibrium point and the 
axes indicated. semi-axes are 4.93 and 


| 
* 
[eo] 
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Fic. 


The Asymptotic Orbit 
The equations for the asymptotic orbit are 


- (i) ~kat 
i=2 k=1 
i 


1.53. 


The initial conditions are chosen that the projection the x-axis the 
initial displacement from the equilibrium point The approximate orbit 
line through the equilibrium point and having the slope 0.92. 


Asymptotic Orbits Near (0.625, 0.625) 


there are purely imaginary roots Indicial Equation (8) for the 
equilibrium point (0.625, 0.625) there are periodic but only asymptotic 
orbits near this point. Since the roots are complex, and not real, the con- 
struction the asymptotic orbits differs somewhat from that used for the 
asymptotic orbits near (0.425, 0.425). 

before begin with Equations with Equations (10’) having 
d/dt. then equate the coefficients the various powers 
Equations and, proper choice the constants integration, 
retain only those terms that are the form where real 
tive and P(t) periodic function there are two undetermined con- 


Hence 


The differential equations the first step are the same the equations 
variation, the general solutions which are 


| 
0.1 
| 
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(26) 


the real part the exponential positive, such terms not 
approach zero approaches infinity and must therefore rejected 
(1) 


putting The desired solutions that satisfy the initial con- 


(27) 
The differential equations this step are 
(D? d,)uy (2D d2)v = 
(2D do) uy (D? = 
(28) 
where are real constants. 
The complementary functions having the asymptotic property are 
where are arbitrary constants. 
The particular integrals, symbolic form, are 
(29) 


where the determinant Equation (8). Now terms explicitly can 
arise Equations (29) only terms occur and and 
such terms appear and obtain the particular integrals 
make use the well known property 


sin sin 
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where f(D) polynomial the operator therefore find the parti- 
cular integrals for and have the same form Then the 
complete solutions that satisfy the asymptotic and initial conditions are 


where are real constants and 
The General Step 


The remaining steps the integration are similar Step terms 
occur the particular integrals. The solutions for and have the 
form 


-(n+Dat | COS cos aes cos 


The expressions the brackets denote sums cosines and sines the 
multiples indicated. The constants integration are chosen that 
The final form the solutions 
n=1 n=1 
consider only the linear terms find approximate orbit 


Table III gives the corresponding values for and for the various values 
listed. 


TABLE III 
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The graph the approximate orbit given Fig. 


© x’? X=. 625 


10. Orbit the Vicinity the Origin 


for the equilibrium point (0.625, 0.625) only one type orbit, viz., the 


asymptotic, exists the vicinity the origin. The solutions for this orbit 


are 


sin 


cos 
sin 


sin 


The higher degree terms are similar those Equations (30). 


Table 


contains the values and for various values for the approximate orbit 


obtained taking only the linear terms Fig. represents the orbit 


thus obtained. 


TABLE 


¥s 625 
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11. Three Dimensional Periodic Orbits 


shall now consider certain three-dimensional periodic orbits the 
vicinity the equilibrium points. 

discuss the motion near libration point make Equations (1) 
the substitutions 

The resulting differential equations are the same Equations (6) with the 
terms added. 

They are 


The equations variation are 


where have the same values Equations (7) and 


The values for the equilibrium points are follows: 


(33) 


(0, (0.465, 0.465) (0.625, 0.625) 


| 
| 
| 
| 
| 
| 
| | 
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The solutions Equations (33) are 


(34) 


where and are constants integration. 
order obtain three-dimensional orbit must take the generating 
solutions either 
period 


(35) 


2 


(0.79 B, 0.61 cos Br 
(—0.61 0.79 By) sin 


where are relatively prime integers. 
the latter case, and must commensurable. this latter case 
the orbits correspond the orbits Class Moulton’s 
doubtful whether and are commensurable, and consequently the con- 
struction these orbits will not considered. 
Expanding Equations (32) power series and cancelling off 
obtain the equations corresponding Equations (10). They are 


4 


higher degree terms 
higher degree terms 


4 
3.2 5 
8 w? 
i=l AS 
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guide the integration Equations (36) power series substitute 


i= i= i= i 


Equations (36) and equate the coefficients the resulting equations, 
cited Equations (36’). 


the first step, differential equations that are the same the equations 
variation are obtained. The periodic generating solutions are 
period (38) 
period 27/y. 
select the orbits having the period 27/8 must put and this 
leads the two-dimensional orbits already constructed. 
that and are commensurable disregarded, and hence the only new orbits 
which will arise are those having the period The solutions the 
first step are therefore 
The desired solutions Step are therefore 
the second step have 


Beginning with Equation (40, 3), the last Equations (40), must have 
Equations (39). The complementary functions Equations (40, and 
must annulled satisfy the periodicity conditions. Hence the solutions 


4 
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this step are the particular integrals Equations (40, and 2). The 
results for the different equilibrium points are listed Table 


TABLE 


(0, (a) (1.255, (b) (0.465, 0.465) (0.625, 0.625) 
1.06 1.52 1.04 
Here the differential equations are 


Table are listed the values for the four equilibrium 


points. 
TABLE 
(0, 01 (a) (1.255, 0) | (b) (0.465, 0.465) (c) (0.625, 0.625) 
| ( - 28: =) sinytT| (— 262 — 2.13) (— 262 + 15.32) (— 262 — 2.39) sin yr 
| | 
128 
histo sin 3yT | + 1.03 sin 3yr | —3.50 sin 3yT — 0.60 sin 3yT 
62 64/27 | —1.06 7.66 —1.20 
— 0.42 sin yr | 0.33 sin yT — 0.57 sin yt — 0.21 sin yr 
+ 0.14 sin 3yr | —0.11 sin 3yT + 0.19 sin 3yr + 0.07 sin 3yr 


The remaining integrations can made the preceding steps. 

Near the origin the and and the odd w’s are all zero. The terms 
consist sums sines odd multiples the highest multiple being 


4 1 
S 
t=1 A i 
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the vicinity the other three libration points considered, the and 
v's with even subscripts and the and with odd subscripts are all zero. 


Libration point (1.255, 


a,2k 


Libration point (0.465, 0.465) 


(41) 


where the C’s and S’s are constants. 


Libration point (0.625, 0.625) 
The solutions are similar (41). 


substituting these results (38) and Equations (31) 


obtain the following: 
Libration (0, 


7=0 k=0 


Libration point (a) (1.255, 


1 


i=0 k= 


Libration point (b) (0.465, 0.465) 


(27 + 2) 9749 
1=0 k=0 


Libration point (0.625, 0.625) 


| 
itl 
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12. Approximate Orbits 


Libration point (0, 


This orbit lies wholly the 


Libration point (a) (1.255, 


1.255 (0.58 0.66 cos 2y7) 


period 


The projections the orbit the planes reference are follows: 


xy-plane 


(0.66 


The major axis 


This ellipse with semi-axes, 0.66 and 0.8 
parallel the y-axis. The graph given Fig. 


The orbit parabola, Fig. 


The orbit double loop, Fig. 


The composite orbit shown Fig. but attempt has been made draw 


the diagram accurate scale. 


The spatial curve figure-8 bent around tangential elliptic cylinder. 
The axis the cylinder vertical and the major axis its elliptical cross- 


section parallel the y-axis. The tangential cylinder encloses Point (a). 


period 

| 
Z 
| 
| } 
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Libration point (b) (0.465, 0.465) 


y x(—7), 


the axes are rotated through 45° the substitutions 


the projections the co-ordinate planes reference are follows: 


The orbit ellipse with the major axis parallel the x’-axis and the 
centre the right 


The orbit parabola similar Fig. 


The orbit double loop similar Fig. The composite result 
figure-8 curve somewhat similar the one shown Fig. 


y 


Fic. 


The figure-8 curve symmetrically situated with respect the diagonal 
the square which the finite masses are the corners. The tangential 


~ 
| 
\ 
| 
| 
, 
Y 
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cylinder vertical and encloses Point (6). The major axis its elliptical 
cross-section along the diagonal the square. The point where the figure-8 
curve crosses the diagonal lies the side (6) remote from the origin. 


Libration point (0.625, 0.625) 


Rotating the axes through 45°, find the projections the planes 


The orbit ellipse lying wholly between the origin and (c). 


This parabola Fig. but with the vertex between the origin and (c). 


The orbit double loop Fig. 

The spatial curve figure-8 the previous cases. The tangential 
cylinder has its major cross-section axis along the diagonal the square 
having the masses corners. The cross-section the tangential cylinder 
formed the xy-plane lies between the origin and Point (c). 


The corresponding spatial orbits near the other libration points can 
readily found rotation the xy-axes through multiples 45°. 


13. Convergence the Solutions 


Only the formal construction the periodic and asymptotic solutions 
has been the solutions are power series parameter, the question 
their convergence arises. 


With respect the periodic solutions, their convergence may established 
theorem (22, Vol. pp. 338-343) showing that periodic solutions exist and 
will converge for all values the time, provided sufficiently small numer- 
ically. Alternatively, may use theorem (19) where 
shown that, the formal construction periodic solutions can made, such 
solutions will converge for all values the time, provided the parameter, 
terms which the solutions are expanded, sufficiently small numerically. 
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The convergence asymptotic solutions was treated Poincaré (22, 
Vol. pp. 338-343). showed that such solutions will converge 
approaches provided the real part the characteristic exponent terms 
which the solutions are expanded real and negative. This condition 
satisfied the case the asymptotic orbits that have been constructed. 
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